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ABSTRACT 

Using rcdshifts as proxy for galaxy distances, estimates of the 2D transverse peculiar velocities 
of distant galaxies {cz ^ 2 x 10^ kms~^) can be obtained from Gala's measurements of proper 
motions. Owing to the large number of galaxies expected to be observed by Gaia, these transverse 
velocities are likely to supersede traditional probes of the large scale velocity field based on 
current and future distance indicator measurements. This Gaia probe of large scale motions is 
completely independent of any intrinsic relations between galaxy properties, hence it is essentially 
free of selection biases. It is also free from homogeneous and inhomogeneous Malmquist biases 
that typically plague distance indicator catalogs. Further, it provides additional information to 
traditional probes which yield line-of-sight peculiar velocities. 

Subject headings: Cosmology: large scale structure of the Universe, dark matter 



1. Introduction 

In the standard cosmological paradigm, pecu- 
liar motions (i.e. deviations from Hubble flow) 
of galaxies are the result of the process of gravita- 
tional instability with overdense regions attracting 
material, and underdense regions repelling mate- 
rial. The coherence and amplitude of galaxy flows 
are a direct indication of the distribution of the 
dark matter, the cosmological background, and 
the underlying theory of gravity Traditionally the 
peculiar velocity field is derived from observations 
of dist ance indicators such a s the TuUy-Fisher re- 
lation (jTuUv &: Fisher! 119771 ) between luminosity 



and rotational velocity of galaxies. The observed 
flux and rotational velocity are then used to infer 
the distance from the TF relation. The distance 
is then subtracted from the redshift, cz, in order 
to obtain the line-of-sight component of the pecu- 
liar velocity of a galaxy, with a typical la error 
~ 0.2 cz. 

Here we point out an alternative probe of the 
large scale velocity field. The Gaia (Global As- 
trometric Interferometer for Astrophysics jl| space 
astrometric mission will likely observe 10^ — 10^ 
galaxies down to V ^ 20 — 25 dur ing the pe- 
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riod 2013-2018(e.g. iRobin et al.l l2012[). Gaia will 
achieve a ~ lO^asyr"^ end-of mission accuracy 
in measurements of proper motions of objects at 
V = 15. Distances to galaxies are needed to derive 
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their transverse peculiar velocities (in kms"'^) 
from the proper motions. However, Gala's par- 
allax measurements give extremely poor distances 
for galaxies at 10s of Mpc away from the observer. 
Therefore, as a proxy for the distance we will use 
the galaxy's redshift, cz, which differs from the ac- 
tual distance by the radial peculiar velocity. The 
relative error in the transverse velocity as a re- 
sult of this approximation is small and decreases 
with redshift. Gala's cnd-of-mission accuracy 
at V — 15, leads to ^ 0.6 cz error in transverse 
velocities. This is significantly larger than the un- 
certainty in line-of-sight peculiar velocities from 
distance indicators. However, we will show that 
the large number of galaxies expected to be ob- 
served with Gaia will beat the increased scatter, 
possibly making Gala's proper motions an excel- 
lent probe of the large scale flows. This probe of 
the large scale flows is completely independent of 
any assumption on the intrinsic relations of galax- 
ies. Further, the 2D transverse motions are or- 
thogonal (in information content as well as in ge- 
ometry) to standard line-of-sight peculiar veloci- 
ties. 

The outline of the paper is as follows. In fj^l 
we present the general set up and describe theo- 
retical tools for analyzing future transverse veloc- 
ity data. We present, in fj3l a rough estimate of 
the expected error in the transverse velocity ob- 
tained by smoothing individual velocities. In the 
concluding section gl we present a general assess- 
ment of the transverse velocity data in compar- 
ison to other probes of large scale motions. We 
also discuss possible sources for redshifts of the 
population of galaxies expected to be observed by 
Gaia. 

2. Methodology 

We will assume an all-sky catalog of redshifts 
and proper motions. We denote the physical pe- 
culiar velocity by v and the real space comoving 
coordinate by r, both expressed in kms~^. Fur- 
ther, v\i = V r and v± =v — v^^r are, respectively, 
the components of V parallel and perpendicular to 
the line-of-sight, where f is a unit vector in the 
line-of-sight direction. We restrict the analysis to 
cz ^ 15000 kms^^ and neglect cosmological geo- 
metric effects, so that the redshift coordinate is 
s = r + v\\r. Note s = r and cz = r + v\\ = s r = s. 



Proper motions transverse to the line-of-sight will 
be denoted by fi. The transverse 2D space velocity 
is of a galaxy at real space distance r is 



r/jb 

677.22 



(1) 

70 kms-^Mpc"^ 



1/iasyr ^ 10^ kms ^ 



which corresponds a transverse peculiar velocity of 
474 kms"^ for l^asyr"^ at d = 100 Mpc. Here- 
after we use Ho = 70 kms~^Mpc~^ to set the 
distance scale. 

However, the distances, r, are unknown, and, 
therefore, we make the approximation 



(2) 



This introduces a relative error uy/cz in the de- 
termination of vj_ where < >i/2^ 200 — 
300 kms^^ (?). Hence the error is negligible as 
we go to cz ^ 2000 kms^^. The error is also ran- 
dom since < v±v\^ >= 0. 

Therefore, the estimated velocity field will be 
given as a function of the redshift space coordi- 
nate. To linear order, velocity fields expressed 
in real and redshift spaces are equivalent. In 
the quasilinear regime, dynamical relations can 
be derived for the velo city field in redshift space 
(|Nusser fc Davislll995[) . ;hanks to the interesting 
property that an irrotational (or potential) fiow 
in real space remain irrotational a lso in redshift 
space ( Chodorowski fc Nussei 1999t ). 



2.1. 



From 2D transverse velocities to 3D 
flows 



Here we offer basic expressions for the deriva- 
tion of the full peculiar velocity field v{s) from 
the smoothed 2D transverse velocity field, v±(s). 
Assuming a potential flow v{s) = — V<E'(s) and 
expanding the angular dependence of $ in sphcr- 
i cal harmonics, ^(s) = $/m(s)>lm(s), gives 
()Arflcen fc Webeijl2005l) 



Im 



ds 



«± = - 2^ — *;,n 



(3) 



(4) 



where '4';„i ~ ^VY/m is the vector spherical har- 
monic. Thanks to the orthogonality conditions 



2 



coefficients can be recovered by 



^)^ii'^mm' tl^c potential 



1(1 + 1) 



(5) 



for I > 0. This means that $(s) can be recovered 
from the v± up-to a monopole term which corre- 
sponds to a purely radial flow with zero transverse 
motions. That is not a serious drawback since the 
monopole term can always be removed from the 
predictions of any model to be compared with the 
data. 

2.2. Testing the potential flow anzats 

Initial conditions in the early Universe might 
have been somewhat chaotic, so that the original 
peculiar velocity field was uncorrelated with the 
mass distribution, or even c ontained vorticity (e.g. 
Christopherson et al.l 120111 ). At late time, a cos- 
mological velocity field should have a negligible ro- 
tational component, 1^'°* on large scale, away from 
orbit mixing regions. The reason is that any cir- 
culation, r = ^ • ds, is conserved by Kelvin's 
theorem. Hence, any rotational component will 
decaying as 1/a, where a is the scale factor. In 
contrast, the irrotational component of the pecu- 
liar velocity will have a growing v ~ -^/a. There- 
fore, on large scales, away from collapsed objects, 
the irrotational component is expected to be neg- 
ligible. The absence of any significant large scale 
vorticity is, therefore, a strong prediction of the 
standard cosmological paradigm. To assess this 
prediction, the observed transverse motions can 
be used to constrain the amplitude of the irrota- 
tional component. This can be d one by writing the 
trans verse component of v'^°*' as (jArfken fc Weber 
20051 ) 



Im 1 



(6) 



where = S x VYim belong to another class of 
vector spherical harmonics that satisfy the same 
orthogonality conditions as '4'. Hence, Vj™' is 
equal to the r.h.s of Eq. [5] but with instead 
of Further, / dfi*;„i • ~ 0, hence the 

recovery of the rotational mode is formally inde- 
pendent of the potential fiow mode. 



3. The expected errors 

We provide estimates of the expected random 
errors in the smoothed transverse velocity field, 
v± (s) , as a function of distance from the observer 

B 

The expected la error, cr^, in the measure- 
ment of an object's proper motion depends on 
its G magnitude and, to a lesser ex t ent, o n its 
color (|de Bruiind[2012l) . Ijordi et al.l (|2010f ) give 
relations between Gala's 0(350 — 1000 nm) mag- 
nitude and the more familiar V and Ic- Using 
these relations we find G ^ V — 0.27, which 
is obtained with V — Ic ~ 1 as is appropriate 



for galaxies (IFukugita et al.l 119951 ). Hereafter, we 
will use (Tu as a function of G, according to 



the expression referenced in Ide Bruiind (j2012[ ) as- 
suming a constant V — Ic = 1 for all galax- 
ies. This cr^ is plotted in left panel of Fig. [T] 
Another useful quantity is the average galaxy 
measurement error as a function of redshift, cz, 
i.e. / dmn(m, c2;)ct^(to)/ / dmn(TO, cz) 
where n(m, cz) is the number density of galax- 
ies (at redshift cz = s) per unit volume and 
per magnitude m and is the magnitude 

cut. Given the luminosity function N{M) we 
have n{m, cz)dm — N{M)dM where M = m — 
51og2o(cz) — 15 is the absolute magnitude. For 
Gala's G band, we approximate N{M) by the 
Schechter form of the V band luminosity func- 
tion taking into account that G — V = —0.26 as 
explained above. He re we adopt the S chechter 
parameters given by (jBrown et al.ll200ll ). Other 
choices for the Sch echter parameters of th e V lu- 
minosity function (jMarchesini et al.ll2012[ ) do not 
change the results significantly at the magnitude 
limits considered here. The results for the aver- 
age error are shown in the right panel for three 
magnitude cuts. The fiatncss of the curves for all 
magnitude, is a refiection of the fact that number 
of galaxies increases strongly with magnitude. 

Given individual measurements v±i ~ v±{ai) 



^ In generating the smoothed Vj_ (s) care must be employed 
since the transverse directions of galaxies in different sight- 
lines within a filtering window do not point in the same 
direction. This difficulty could b e overcome by ten sor win- 
dow smoothing a la POTENT jPekel et al.Hl990h . How- 
ever, we will not be concerned with these fine details at 
this stage. 
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write the smoothed velocity as 



(7) 



where the summation is over all galaxies, a±i = 
safj,i and is a smoothing window function. 
The 1(7 errors on is given by 



(8) 



The summation over galaxies can be transformed 
into a volume integration with the same argument 
but multiplied by the number density of galaxies. 
Doing so for a uniform distribution and assuming 
a Gaussian window, W ^ of width Rq, we get 



87r3/2i?3, /™'"" dmn(m, cz)ajl'^{m) 



(9) 

we have assumed the distant observer limit so 
that |si — s| ^ s. For a Top-Hat window of 
the same width we get the same expression but 
with 47r/3 as the numerical factor in the denomi- 
nator. Substituting a^im) (see left panel Fig.[T]) in 
Eq. [9l we compute the expected error, ux, in the 
smoothed Vi_ , for a gaussian smoothing with Rq = 
1500 kms"^. The top panel in Fig.[2]shows curves 
of (Jx as a function of distance for three magnitude 
cuts. For comparison the figure also plots the er- 
ror in the filtered line-of-sight peculiar velocities in 
the SFI-I- + catalog of TF mea surements of 4000 
galaxies ( Masters et al. 2006() . There is a signifi- 
cance decrease in tT_L as the magnitude is increased 
from G = 14 to 15, but the improvement is not as 
dramatic when fainter galaxies with 15 < G < 16 
arc included. The reason is the rapid deteriora- 
tion in cr^ at G = 16 which is not compensated by 
the added number of fainter galaxies. At redshifts 
cz ^ 6000 knis"^ and for G < 15, peculiar veloc- 
ities from Gala's proper motions are expected to 
fair much better than the SFI-|--|- catalog. An- 
other important quantity which can be computed 
from transverse velocities is the dipole motion (i.e. 
bulk flow) of spherical shells of a given thickness. 
This motion is described by a constant term B 
and gives rise to a transverse velocity field of the 
form v±_B = B — f(B • f). The dipole term B 
can be found by least squares fitting of v_lb to 
the observed velocities v^_i . The expected error in 



B as a function of distance of the shell, is plot- 
ted in the bottom panel in Fig. [5] Predictions for 
3 magnitude cuts are plotted for spherical shells 
of 3000 kms^^ in thickness. For c omparison we 



also plot the WMAP7 ACDM model (Larson et al 



l201^ predictions for the amplitude of the veloc- 
ity dipole on spherical shells. It is encouraging 
that the predicted amplitude is larger than the 
expected error out to relatively large distances. 

4. Discussion 

The number of galaxies expected to be ob- 
served by Gaia is likely to exceed standard dis- 
tance indicator data by t wo orders of magni- 
tude (iMasters et all l2006t ISpringob et all I2012 : 
Masters et al.l 120081: [Courtois et al.l l201l[) . De- 
spite the larger object-by-object error, this makes 
the method presented here a serious contender to 
traditional probes of the peculiar velocity field. 
The method has several advantages. Firstly, it is 
completely independent of any assumed intrinsic 
relations of galaxies and, hence, it does not suffer 
from the usual concerns related to these relations, 
e.g. linearity, selection biases and dependence on 
environment. Secondly, it yields the 2D transverse 
velocity component and hence it offers completely 
orthogonal information to standard probes which 
yield the line-sight-component. Thirdly, it is free 
from h omogeneous and inhomo geneous Malmquis 
biases (|Lvnden-Bell et al.lll988l) . 

The usefulness of the method for probing the 
3D velocity field on scales of a few 10s of Mpc is 
limited to cz ^ 10^ kms~^. However, large scale 
moments of the velocity field can be assessed at 
much larger distances. In particular the error of 
the dipole on (i.e. bulk flow of) spherical shells 
can be estimated with ~ 100 — 200 kms~^ error at 
cz ^ 1.5 X 10^ kms^^. At larger redshifts, neither 
this method nor traditional ones are comparable to 
the constraints on the dipole from gal axy luminosi 



ties i n future galaxy redshift survevs (|Nusser et al 



20111) . At lower distances (< 2000 kms'^) the the 
transverse motions of galaxies could play an im- 
portant role at providing new constraints on the 
motion of Local Group of galaxies. 

Gaia will provide spectroscopic information of 
unresolved galaxies (?). But the inferred redshifts 
may not be sufficiently accurate for a ll unresolved 
galaxies with astrometric data (e.g. iRobin et al 
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Fig. 1. — Expected error in Gaia's proper motion measurements. Left: la error of an object as a function of 
its G magnitude. Right: mean la error of objects as a function of distant for three values of the G magnitude 
cut. 



20121). The Two M ass Redshift Survey (2MRS; 
Huchra et all l201lh offers redshifts of - 4 x 10^ 
galaxies down to Kg = 11.75. This is the deep- 
est all-sky redshift catalog currently available. It 
was originally planned to reach up to K=12.2 
mag and to include ~ 10^ galaxies. Using G « 
V — .27 and the observe d V band luminosity func- 
tion (jBrown et al.l 120011 ) . the expected number of 
galaxies brighter than V = 15.27 (i.e. G ~ 15) 
is ~ 10^, similar to the expe cted number for 
Ks = 12.2. Further, Ks ^ K dCarpenteilbOOll) 
and V — K > 2.7 for most galaxies ( AaronsonI 
1978h . we conclude that the = 12.2 2MRS red- 



shifts would correspond to of most of Gaia galaxies 
observed to G = 15. This is particularly inter- 
esting as Gaia's astrometric accuracy deteriorate 
rapidly at fainter objects, however However, it is 
unclear if 2MRS will be continued to Kg = 12.2 
in the very near future (Macri, private commu- 
nication). For the purpose of the analysis pre- 
sented here one could use a catalog of photomet- 
ric redshifts based on the 2MASS galaxy catalog 
(Skrutskie et al. 2006), containing almost 1 mil- 
lion sources with Kg < 13.5 mag. Its current 
form (2MASS XSCz, Jarrett 2004) offers errors as 
large as 20-25%, which will improve in the coming 
years using the data from other galaxy catalogs for 
the photo-z estimation (Bilicki, private communi- 
cation). 



We have restricted the error analysis here to 
G '--^ 15 since redshifts will probably not be avail- 
able for all fainter galaxies . However, data at 
fainter magnitudes can well be exploited by com- 
puting the dipole clS Bi function of an effective 
depth corresponding to a certain magnitude range. 
This can then be compared with model predict ions 
for an equivalent quantity (jBilicki et al.ll201ll ). 
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